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Abstract 

The air-fluorescence yield is a key parameter for determining the 
energy scale of ultra-high-energy cosmic rays detected by fluorescence 
telescopes. A compilation of the available measurements of the abso- 
lute air-fluorescence yield normalized to its value in photons per MeV 
for the 337 nm band at given pressure and temperature has been re- 
cently presented in |1J. Also, in that paper, some corrections in the 
evaluation of the energy deposited in the corresponding experimental 
collision chambers have been proposed. In this note this compari- 
son is updated. In addition, a simple statistical analysis is carried 
out showing that our corrections favor the compatibility among the 
various experiments. As a result, an average value of 5.45 ph/MeV 
for the fluorescence yield of the 337 nm band (20.1 ph/MeV for the 
spectral interval 300 — 420 nm) at 1013 hPa and 293 K with an un- 
certainty of 5% is found. This result is fully compatible with that 
recently presented by the AIRFLY collaboration (still preliminary) in 
such a way that including this latest result could even lowered the 
final uncertainty below the 5% level with high reliability. 



The fluorescence technique has been proved to be very fruitful for the 
study of ultra-high-energy cosmic rays. The air-fluorescence yield, defined as 
the number of fluorescence photons per unit deposited energy, is a key cali- 
bration parameter of the fluorescence telescopes. In we presented a com- 
parison between the laboratory measurements of the absolute air-fluorescence 
yield available in the literature normalized to common units (photons/Me V), 
air conditions (dry air, 800 hPa and 293 K) and wavelength interval (337 nm 
band). From the comparison of this normalized 1337 parameter, discrepan- 
cies larger than the quoted uncertainties of the measurements were found. 
According to the Monte Carlo analysis carried out in pQ, non- negligible cor- 
rections in the evaluation of the energy deposition should be applied to these 
results, in particular in experiments where authors neglected the energy de- 
posited outside the field of view of the light detector. In fact, the proposed 
corrections are even larger than the quoted uncertainties in some cases (i.e., 
measurements of Kakimoto et al. [2] and Lefeuvre et al. [3]). As a result of 
these corrections, the Y337 values get closer to each other in such a way that 
they are contained within the 6 — 7 ph/MeV interval at the above reference 
pressure and temperature [1]. 

In this work, the calculations presented in [1] have been updated and the 
compatibility of the fluorescence yield results with and without our correc- 
tions has been evaluated quantitatively. This has allowed to obtain a reliable 
1^337 value resulting from an average of the available absolute measurements. 
With respect to the MC analysis performed in [1] , some improvements in the 
algorithm for the evaluation of the energy deposited in the collision chamber 
of the experiments have been made. In particular, the density correction 
on the cross section for the ionization processes leading to ejection of K- 
shell electrons, which was neglected in |T], have been included in this work. 
These highly excited ions generate X rays which deposit all their energy 
within a very short distance. As a consequence, the density correction for 
the K shell lowers the energy deposited by high-energy electrons (GeV range). 
GEANT4 simulations of energy deposition in air carried out by MACFLY |1] 
and AIRFLY [S] are in agreement with our calculations at the level of 2% 
(1%) for electron energies in the GeV (MeV) range. 

In table [l| we present the compilation of fluorescence yield values [1] 
with these updated correction^ In this case the fluorescence yield have 

^The updated corrections are very similar to those proposed in [1] except for experi- 
ments working with very-high-energy electrons; however, the average fluorescence yields 
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Table 1: Comparison of available measurements of the absolute air- fluorescence yield 
normalized to the 337 nm band at 1013 hPa and 293 K (dry air). Experiments are listed in 
column 1, the electron energies are shown in column 2 and the corresponding results, both 
uncorrected and corrected (in bold), are listed in column 3. The experimental uncertainties 
quoted by the authors are shown in column 4 and the size of the proposed correction is 
displayed in the last column. 



Experiment 


E (MeV) 


Ys37 (ph/MeV) 


Quoted error 


Correction 




1.4 


4.5 / 4.8 




+6% 


Kakimoto [2j 


300 
650 


4.4 / 5.5 
3.8 / 4.8 


10% 


+25% 
+27% 




1000 


4.3 / 5.5 




+29% 


Nagano |£j 


0.85 


5.0 / 5.4 


13% 


+6% 


Lefeuvre j3j 


1.1 


5.1 / 5.5 


5% 


+7% 


1.5 


5.6 / 6.1 




+8% 




1.5 


4.3 / 4.4 




+1% 


MACFLY [4J 


20 • 10^ 


4.4 / 4.3 


13% 


-2% 




50 • 10^ 


4.6 / 4.5 




-2% 


FLASH [8J 


28.5 • 10^ 


5.5 / 5.6 


7.5% 


+2% 


AirLight |9] 


0.2 - 2 


5.8 / 5.4 


16% 


-7% 


AIRFLY [6J 


120 • 10^ 


5.6 / - 


<5% 





been normalized to its value at 1013 hPa (instead of 800 hPa as presented 
in d]) and 293 K. The latest result of AIRFLY [U] carried out using 120 GeV 
protons have been included in this compilation, but the impact on the average 
has been studied separately, as this value is still prehminary. 

The fluorescence yield results have been displayed as a function of en- 
ergy in figures [l] (uncorrected values) and [2] (corrected ones). As can be 
appreciated, they are in a better agreement when including our corrections. 
In addition, the corrected results give more support to the expected energy 
independence of the fluorescence yield. 

In order to know quantitatively to what extent our corrections favor the 
agreement between the absolute results included in this comparison, a statis- 

presented below are basically not affected by this update. 
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Figure 1: Comparison of uncorrected normalized 1337 values as a function of energy. The 
horizontal line represents the corresponding weighted average value (see text for details). 




♦ Kakimoto 

♦ Lefeuvre 

♦ FLASH 

— <Y>=5.48 ph/MeV 



♦ Nagano 
MACFLY 

♦ AirLight 

-•-Th. prediction (scaled) 



100 

E (MeV) 



1000 



10000 



100000 



Figure 2: Same as figure [l] for the corrected 1337 values. The black circles connected by 
lines represent the weak energy dependence predicted by our simulation |10j . 
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Table 2: Comparison of absolute Y^^t values of table [ij where measurements reported 
by a given experiment at different electron energies have been averaged. 



Experiment 


^3™' (pli/MeV) 


Y,%r (ph/MeV) 


Quoted error 


Kakimoto |2j 


4.3 


5.2 


10% 


Nagano [7] 


5.0 


5.4 


13% 


Lefeuvre [3] 


5.4 


5.8 


5% 


MACFLY [4J 


4.5 


4.4 


13% 


FLASH [8] 


5.5 


5.6 


7.5% 


AirLight [9] 


5.8 


5.4 


16% 


AIRFLY [g 


5.6 




<5% 



tical analysis has been performed. In the first place, for a given experiment, 
results obtained at different energies have been averaged assuming a common 
systematic erroi|^ The results are shown in table [2] for both the uncorrected 
(column 2) and corrected (column 3) 1^337 values together with the percentage 
uncertainty associated to each experiment (column 4). 

In the second place, the weighted average of the data sample shown in 
table |2] has been computed from 

{Y) = , (1) 

using as weights the reciprocal of the quoted square uncertainties (i.e., Wi = 
I /erf). If these af are assumed to actually represent the variances of the 
corresponding (normal) probability distributions, our weighted mean would 
be the best estimator of the Y337 value with a variance given by ^/'^i) ^ ■ 
However, note that authors usually do not include any error contribution 
from the evaluation of the deposited energy, and thus, uncertainties are very 
likely underestimated, at least in some experiments. In fact, the statistic 
divided by the number of degrees of freedom (x^/ndf) is found to be some- 
what larger than 1 for the uncorrected sample. Therefore, following the usual 



^This way, no higher significance is given to experiments measuring at several energies 
than those measuring at a single energy. 
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procedure, the variance of the weighted mean was corrected multiplying by 
X^/ndf, that is, 



Data listed in table [2] for both the uncorrected and the corrected sample 
are plotted in figure |3| The average value for the uncorrected Y^^j sample 
(figure g is found to be {Y) = 5.12 ph/MeV with a(y) = 0.23 ph/MeV 
(4.4%) and a x^/ndf value of 1.60, whereas the corresponding results for the 
corrected sample (figure 3b) are (y)=5.48 ph/MeV, a(Y) = 0.20 ph/MeV 
(3.6%) and x^/ndf=1.07. Therefore, our corrections lead to a more consis- 
tent data sample suggesting that they do improve the determination of the 
deposited energy for the different experiments. 

In the absence of other objective criteria, data have been weighted using 
the quoted experimental uncertainties, as mentioned above. Nevertheless, 
this weighting procedure has no significant effect on the final result. For 
instance, assuming the same weight (i.e., the same variance) for all the ex- 
periments the results are (Y) = 5.09 ph/MeV and o"(y) = 0.25 ph/MeV 
(5.0% error) for the uncorrected data sample and (Y) = 5.30 ph/MeV and 
o"(y) = 0.20 ph/MeV (3.8% error) for the corrected one. The weighted sample 
mean turns out to be somewhat larger than the non-weighted mean mainly 
due to the high statistical significance of the results of Lefeuvre et al. and 
FLASH, with the smallest associated uncertainties. In any case, a similar 
decrease in the standard deviation is found when corrections are applied 
whichever weights are chosen. 

As a further test, it has been checked that the weighted mean of the 
corrected data sample does not change significantly if some measurement 
is excluded. For instance, removing the MACFLY result, which is the one 
showing the largest departure from the average, would lead to 5.61 ph/MeV. 
Instead, excluding the result of Lefeuvre et al., with a claimed uncertainty of 
5% that is incompatible with our correction of 7%, an average fluorescence 
yield of 5.23 ph/MeV is obtained. 

We have also studied the energy (in) dependence of the fluorescence yield 
in this data sample. Several experiments have supported the assumed inde- 
pendence of the fluorescence yield at the level of about 10% in several energy 
intervals. On the other hand, the theoretical analysis described in [10] pre- 
dicts a slight increase (~ 2%) with decreasing energies in the 0.1 — 10 MeV 
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(a) Uncorrected Y337 values. 
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(b) Corrected ¥337 values. 

Figure 3: Graphical representation of the ¥337 values at 1013 hPa and 293 K (see 
table [2]). The weighted mean value {Y) (vertical continuous line), the standard deviation 
of the mean a(^Y) (half of interval between dashed lines) and the statistic normalized 
by the number of degrees of freedom are shown in the legends, (a) Discrepancies between 
the uncorrected values are sometimes larger than quoted errors, resulting in a relatively 
large x^/ndf value, (b) Better agreement and consistency of the data sample are found 
when applying the proposed corrections. 
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Figure 4: The same as figure 3b but previously scaling to 100 MeV using the weak energy 
dependence predicted by our simulation (see figure [2]). The consistency of the data sample 
is strengthened slightly further. 



range, which is also compatible with experimental data (see figure [2]). In 
principle, this effect should be included in order to quantify the consistency 
of the available results. Average values listed in table [2] have been recal- 
culated by previously scaling all the measurements to a common electron 
energy of 100 MeV according to the energy dependence predicted by our 
simulation. The results are shown in figure |4j This energy scaling slightly 
lowers the x^/ndf value (from 1.07 down to 1.00) while both the average 
value and its uncertainty remain nearly unchanged (from 5.48 ± 0.20 to 
5.43±0.19 ph/MeV). Although this result would support the weak energy de- 
pendence predicted by our simulation, we note that the evidence is still very 
weak and this small reduction of x^/ndf after the energy correction might 
have happened by chance. Therefore, even if confirmed experimentalljj^ this 
weak energy dependence would not have any relevant impact on the energy 
reconstruction of cosmic rays using the fiuorescence technique. 

Taking into account all the above considerations, including the possible ef- 
fects of the weak energy dependence of the fiuorescence yield and the weight- 
ing procedure, a conservative result from our analysis is Y337 = 5.45 ph/MeV 



•^Currently not possible due to the accuracy limitations in the determination of the 
energy dependence of the fluorescence yield in this large range. 
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with an estimated error of 5%. According to the comparison of our simula- 
tion of the energy deposition with GEANT4, a small systematic uncertainty 
of below 2% should be added, although it does not affect the found improved 
compatibility of results (i.e., the values). The recent absolute measure- 
ment of the AIRFLY collaboration [H] yields Y337 = 5.6 ph/MeV with an 
uncertainty of < 5% (still preliminary), which is fully compatible with the 
above value. If this new result is included in the average, then a weighted 
mean value of 5.52 ph/MeV is obtained from ([T| with an uncertainty of < 5%. 

As already mentioned, for the comparison presented here we have nor- 
malized the air-fluorescence yield measurements to its value for the 337 nm 
band at 1013 hPa and 293 K. However, in some occasions it might be more 
convenient to use the integral of the fluorescence yield in a wider spectral 
range and/or other pressure and temperature conditions. The conversion 
can be easily done following the procedure described in detail in [H [10]. For 
instance, the above average value would be of 20.1 ph/MeV (ib5%) for the 
300 — 420 nm spectral range at the same reference pressure and tempera- 
ture, which would become 20.3 ph/MeV if the measurement of AIRFLY is 
included. 

Our simulation can also provide a theoretical value of the air-fiuorescence 
yield. Unfortunately, the evaluation of the fluorescence emission cannot be 
very precise due to the large uncertainties in the relevant molecular param- 
eters. Therefore, we expect a large uncertainty in this calculation of the 
fluorescence yield, which we estimated to be about 25% [I^. Nevertheless, a 
result]^ for I337 of 6.3 ph/MeV has been found, which is consistent with the 
experimental ones, providing a valuable theoretical support to these mea- 
surements. 

In summary, the corrections to the measurements of the absolute air- 
fluorescence yield discussed in p!] increase significantly the compatibility of 
the various experimental results. An average value of 1337 = 5.45 ph/MeV 
with a 5% uncertainty has been obtained. If the absolute fluorescence yield 
and error of AIRFLY are confirmed, a consensus on this important parameter 
with an uncertainty below the 5% level could be reached with high reliability. 

''Using the quenching parameter of pT| . 



8 



Acknowledgements 

This work has been supported by the Spanish Ministerio de Ciencia e Innova- 
cion (FPA2009-07772 and CONSOLIDER CPAN CSD2007-42) and "Comu- 
nidad de Madrid" (Ref.: 910600). J. Rosado acknowledges a PhD grant from 
"Universidad Complutense de Madrid". The authors thank our colleagues 
of the Auger collaboration, in particular D. Garcia Pinto and V. Verzi, for 
fruitful discussions and comments on this work. 

References 

[I] J. Rosado, F. Blanco and F. Arqueros, Astropart. Phys. 34 (2010) 164. 
[2] F. Kakimoto et a/., Nud. Instr. and Meth. A 372 (1996) 527. 

[3] G. Lefeuvre et al, Nud. Inst, and Meth. A 578 (2007) 78. 

[4] P. Cohn et al. [MACFLY Collaboration], Astropart. Phys. 27 (2007) 317. 

[5] AIRFLY collaboration, private communication. 

[6] M. Ave et al. [AIRFLY Collaboration], Proceedings of Cosmic Ray Inter- 
national Seminars (Catania, Italy, 2010). Available in arXiv:1101.3799. 

[7] M. Nagano et al, Astropart. Phys. 22 (2004) 235. 

[8] R. Abbasi et al. [FLASH Collaboration], Astropart. Phys. 29 (2008) 77. 
[9] T. Waldenmaier et al, Astropart. Phys. 29 (2008) 205. 
[10] F. Arqueros et al, New J. Phys. 11 (2009) 065011. 

[II] M. Ave et al [AIRFLY Collaboration], Astropart. Phys. 28 (2007) 41. 



9 



